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We study the pairwise concurrences, a measure of entanglement, of the ground states for the 
frustrated Heisenberg ring to explore the relation between entanglement and quantum phase tran- 
sition associated with the momentum jump. The groundstate concurrences between any two sites 
are obtained analytically and numerically. It shows that the summation of all possible pairwise 
concurrences is an appropriate candidate to depict the phase transition. We also investigate the role 
that the momentum takes in the jump of concurrence at the critical points. We find that an abrupt 
momentum change rusults in the maximal concurrence difference of two degenerate ground states. 
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I. INTRODUCTION 

A quantum system changes its groundstate properties 
in a fundamental way when quantum phase transitions 
(QPTs) occur at absolute zero Q, which are induced by 
the change of external parameters or coupling constants 
and are driven by quantum fluctuations. In a finite sys- 
tem the discontinuity of the groundstate energy is often 
used to characterize the occurrence of QPTs for there 
is an energy-level crossing at the critical point. We will 
show, in this paper, that such kind of phase transitions 
are usually accompanied by the change of symmetry char- 
acterized by conservative quantities such as the momen- 
tum which is the generator of translation. Our investiga- 
tions will relate to the quantum entanglement measured 
by the total concurrences 0- 

Actually, quantum critical points are governed by a 
diverging correlation length and there exists a very close 
relation between quantum correlation and quantum en- 
tanglement, which is known as the resource that en- 
ables quantum computation and communication. So ex- 
ploring the role of entanglement in a phase transition 
has attracted great attention from both the commu- 
nities of quantum computation and quantum statistics 
H II H IE II 11 S H O IH 13 03 People have con- 
nected the theory of critical phenomena with quantum 
information by exploring the entangling resources of a 
system close to its quantum critical point. They demon- 
strate, for a class of magnetic systems or the interacting 
quantum lattice spin systems at zero temperature, that 
entanglement shows scaling behavior in the vicinity of 
the transition point where the level crossing occurs at 
degenerate groundstates [lfij - 

From these studies we observe that the degeneracy of 
groundstates at the critical point may result in the uncer- 
tainty of entanglement. This fact means the discontinu- 
ity of concurrence in the vicinity of the phase transition 
point. Such sudden change of the concurrence as the 
variation of an external parameter or coupling constant 



is vividly called the jump of concurrence. However, the 
groundstate energy-level crossing may not always result 
in a jump of a certain type of concurrence, such as the 
NNN or other pairwise concurrences. On the other hand, 
this kind of QPTs must be accompanied by a change of 
a certain conservative quantity, such as the momentum 
or the macroscopic magnetization. This observation may 
provide us a new way to find an appropriate definition 
of concurrence that just characterize the property of the 
quantum spin systems at the critical point. 

In this paper, we consider a frustrated Heisenberg ring 
system which contains rich phases in the ground state. To 
reveal the connection between the concurrence behavior 
and the symmetries of the separated phases around the 
critical point, we study the difference between the con- 
currences of two degenerate ground states at the critical 
points. We find that there does exist such a discontinuity 
of a conservative physical quantity-momentum (the gen- 
erator of translation) which results in the maximization 
of entanglement difference of the two degenerate ground 
states. 

The paper is organized as follows. In Sec. II, numer- 
ical result for the concurrence behavior of a frustrated 
Heisenberg ring around the critical points is given. It 
shows that a single type of concurrence is not sufficient 
to depict a QPT, while the summation of all types of 
concurrence may be. In Sec. Ill, the exact results are 
employed to explain our observation concluded in Sec. II. 
In Sec. IV, the role that another conservative quantity 
- momentum plays in the concurrence jump is studied. 
The summary and some discussions are given in Sec. V. 



II. CONCURRENCE JUMPS IN QPTS 

We start from a ID frustrated spin-i Heisenberg model 
with periodic boundary conditions, which belongs to the 
Majumdar-Ghosh (MG) families of models (see Fig. 
1)- 
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FIG. 1: The schematic structure of a frustrated spin-i 
Heisenberg ring. Jo denotes the nearest neighbor coupling 
constant, while 3\ and J2 denote the next nearest neighbor 
coupling constants between even and odd number of sites re- 
spectively. In this paper we only concern the simple case 
Ji = J 2 = J. 



The Hamiltonian reads 



JY 



N 



H — Jo ^ Si • Si+i + J ^ Si ■ S 



i+2, 



(1) 



where N is even denoting the total number of the lattice 
sites, Si is the spin operator at ith site, and Jo (J) is the 
strength of the NN (NNN) exchange interaction. For the 
periodic boundary conditions, we have Sjv+i = Si. For 
arbitrary Jo and J, the groundstate energy of Hamilto- 
nian (1) can be formally written as the function of the 
parameters Jq and J 



E g (J , J) = J h + Jh, 
in terms of the NN correlation function 

h o = Yl (0I Sl ■ Si + 1 1-9) 

i 

and the NNN correlation function 

h = ^2(g\ S s : • S j:+2 |s) 



(2) 



(3) 



(4) 



with respect to the ground state \g). 

Now we consider the behavior of groundstate energy 
as the variation of parameters J and Jo- For an infi- 
nite system, the energy-level crossing should induce the 
discontinuity of the derivative of groundstate energy 



dE g , , dH 

-jry- = .9 Try .9 = h , 

oJ dJ 



(5) 



which is just the NN correlation function. Since there 
exist an algebraic relationship between pairwise concur- 
rence and correlation function |T^ | 



C H = 77 m ax{0, \G X 



Gyy I G z 



-1}: 



(6) 



where G aa = (g\ af cr" \g) (a = x,y,z) are correlation 
functions, then according the Eq. © the energy-level 
crossing will lead to the discontinuity of the NN con- 
currence. On the other hand, one can also establish a 



similar relation between NNN correlation function h and 
the discontinuity of dE g /8J, 



dE g ,,dH 



(7) 



Notice that, for nonzero Jo, if the correlation h has a 
jump due to the discontinuity of OE g /dJ, the other cor- 
relation function ho must experience a jump at the same 
point J — J c . Actually, at the energy-level crossing point 
the ground states are degenerates, i.e. 



where 



JqIiq + Jh = Joh' + Jh', 

i 



(8) 



(9) 



are the corresponding correlations with respect to an- 
other ground state \g'). No doubt, both the NN and 
NNN correlation functions must be discontinuous or have 
jumps at critical point J = J c . But the jump of one type 
of correlation may not necessarily induce the jump of its 
corresponding concurrence. A natural question is that, 
which type of concurrence, NN or NNN, plays a major 
role in depicting the QPTs at the critical points? To an- 
swer this question, we investigate the frustrated spin-i 
Heisenberg ring numerically and analytically. In Fig. 2, 
the eigen energies of ground and first excited state are 
plotted for the systems with size N = 6, 8, 10 and 12. A 
and B denote the two energy-level crossing points. 

The reduced density matrix for two spins located at 
sites i and j has the form 



/ va 



Wi 



V Vij ) 



\ 






(10) 



with respect to the standard basis vectors Iff), |f \ ), | \ f ), 
and Iff). Correspondingly the concurrence can be calcu- 
lated by 



max{0,2(|^ y | -««)}• 



Furthermore, for the groundstate with vanishing total 
spin, it can be connected to the isotropic correlation func- 
tion (<Tj • <7j) by [l^ | 



Gi- 



max {0, — {pi ■ Oj) — 1} . 



(11) 



Analytical and numerical results show that the ground 
states of the Hamiltonian (1) for the sizes we concern 
have spin zero. Then the concurrences can be obtained 
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FIG. 2: 27ie ground and first excited eigen energies for the 
systems of the size N = 6, 8, 10 and 12. A and B denote the 
energy level crossing points where the quantum phase transi- 
tions occur. The critical point A is always at J — J c = Jo/2, 
while the position of point B depends on the size of system. 

directly from the corresponding correlation functions. 
Now we define the a type concurrence 

JV-l 

CH = Ci, i+ a, (12) 

i 

where a = 1, 2, N/2, and the total concurrence as the 
summation of all types concurrences, i.e. 

N 
2 

c T = J2 cla] - ( 13 ) 

a 

It is obvious that CV shares the same property as the 
average concurrence 19]. We also define the concurrence 
jump as AH = |c[ a] - C [ ° ] \ (A T = \C TL - C TR \) de- 
noting the a (total) type concurrence difference between 
the two ground states at left (L : J = J c — + ) and right 
(R : J = J c + + ) side of the critical points. 

The conjectures of the relationship between entangle- 
ment and QPTs in Refs Q, Q] tell us that the concur- 
rence, a measure of entanglement, should be changed 
largely at the critical points. So we calculate concur- 
rence at point B, for N = 8, 10 and 12 systems numer- 
ically. The corresponding concurrence jumps and 
At are listed in Table. 1. 

The big concurrence jumps AW for N = 8, 10, 12 and 
Al 2 l for N = 8 match the energy-level crossing in Fig. 
1, while the rest concurrences show no special behavior 
at point B. It indicates that, at the critical point, not 
all types of concurrence change largely. It shows that 
different types of concurrences show different behaviors 
around the critical points. It seems that there exists no 



N 8 10 12 

a = 1 0.7660 1.2755 0.6228 

a = 2 1.8228 
a > 2 

total 1.0568 1.2755 0.6288 



TABLE I: The jumps of the pairwise concurrences C^ a ' and 
total concurrence Ct at the critical point B for the N — 8, 10 
and 12 systems. It implies that not all types of, but total 
concurrence are appropriate to depict the QPTs. 

single preferable type of concurrence in characterizing the 
QPTs. A natural way to treat this problem is to use the 
summation of all types of concurrence. In the following 
section, we will investigate the critical behavior at point 
A. The results further demonstrate that the total con- 
currence Ct seems to be a good candidate to depict the 
QPTs. 



III. MOMENTUM JUMP IN QPTS 
CHARACTERIZED BY ENTANGLEMENT 

According to quantum mechanics, we can always find 
out a conservative quantity to distinguish the two degen- 
erate ground states at the energy-level crossing point. In 
this section we study the critical behaviors at the energy- 
level crossing point A and B (Fig. 2) . We will show that 
these critical points are just between the phases with mo- 
menta and 7r. 

In the case J = Jo/2, the exact ground states of the 
N site system can be explicitly expressed as 

\<f>i) = [12] [34] .. [iV — 1 N] (14) 

or 

102) = [23] [45] .. [Nl] (15) 

which are the direct products of the resonant valence 
bond (RVB) states [ij] = (|U) - |1H)/V2 of two spins 
located at the lattice sites i and j |20J. Obversely, \4>i) 
and |02) are not orthogonal except in the case of thermo- 
dynamical limit, but their combinations \<f>i) — \<f>2) and 
\<pi) + 102) are two orthogonal degenerate ground states. 

Because the Hamiltonian is invariant under the trans- 
lational transformation T, where T|T)j = lt),;+ii the 
common eigenstates \vj) = ^(Si, S2,..., Sjy)) of H and 
T has momentum 

27rn 

k = ——- = na, n = 1, 2, N 

which satisfies T\ijj) = exp(ina) \ip). Actually as men- 
tioned above, at point J — Jo/2, one can construct the 
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two degenerate ground states as 

1 



|-0x> = 

IV>2> = 



1 

Mi 



|01> - |02>] , 



(16) 



with momentum fc = and 7r respectively. Here 



«i = 
« 2 = 



(<f> 2 \<j> 2 ) - 2Re{{<t>i |&», 
(02 |^2> + 2 J Re«<^i |^ 2 », 



are the normalization factors. 

We take a system of small size as an analytical illus- 
tration for the momentum jump in QPT. For the system 
of N = 6, the ground states of the Hamiltonian can be 
obtained exactly in the whole range of J as 



IVO ( J > Jo/2) 

(KM + qliMwn; (j < Jo/2) 



(17) 



corresponding to the eigenvalues 



E x = -3(J + J)/2, 

£ 2 = fa - 5/2) J + (1/2 - rj) J, 



(18) 



respectively. Here 



and 



fi.3 
JJ = 



8r] 2 
J - 



- 8r, + 20, 

3 Jo + x/9J 2 - 18JJ + 13J 2 



2 (J - Jo) 



/-^ rp _j_ /jiz /^o _|_ ^4 /jio 



(19) 



2V2 

iTmu>-(i-r)iuuu»], 



(20) 



is the excited state of the Hamiltonian (1) for J = Jo 
with eigen energy and momentum n. It is easy to find 
that there is only one critical point at which the energy- 
level crossing occurs. But for N > 6, there is one more 
energy-level crossing points at J > Jo /2 as illustrated in 
Fig. 2. 

Based on the exact results, the two spin concurrences 
are obtained as 



-4 ( V 2 + 2 V - 2) /0 3 (J < Jo/2) 
(J > Jo/2) , 



(21) 



and 



(J < Jo/2) 
(J > Jo/2) 



which are plotted in Fig. 3 



(a = 2,3), (22) 
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FIG. 3: Various types of concurrence as the function of J for 
N = 6. The momenta of the ground state for J < Jo/2 and 
J > Jo/2 are ir and respectively. It shows that only NN 
type of concurrence possesses a jump around the critical point 
J = Jo/2. 



It shows that the NN concurrence has a jump at 
the critical point A, while other types of concurrences 
C^ a \ (a > 1) do not have any special behavior. Obvi- 
ously Ct also has a jump at the critical point, which is 
in agreement with the observation presented in Sec. II. 
In the following section, we will investigate the relations 
between various types of concurrence jumps and a con- 
servative quantity, momentum which has a jump for the 
QPTs. 



IV. CONCURRENCE JUMP IN ASSOCIATION 
WITH MOMENTUM JUMP IN QPT 

From the above analysis, we know that the jumps of 
concurrences C' Q ' at the critical points may be induced 
by the energy-level crossing or the discontinuity of the 
groundstate energy as a function of the coupling con- 
stants. On the other hand, at the energy-level crossing 
point, the ground states are degenerate. Thus an arbi- 
trary linear combination of two degenerate ground states 
is also the ground state. If a certain type of concur- 
rence has a jump at the critical point, the corresponding 
concurrence of the combined ground state should be un- 
certain. Meanwhile the difference of the concurrences be- 
tween the two orthogonal combined ground states should 
also depend on the way of the combination. On the other 
hand, as the energy-level crossing there must exist a con- 
servative quantity which also experience a jump. Then 
the phase separation can be also well described in asso- 
ciation with the jump of such a quantity. In this paper, 
this conservative quantity is momentum which is the gen- 
erator of translation. In general, one may say that the 
discontinuity of 3E g /dJ leads to the jump of concur- 
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rence at the critical points, but on the other hand, one 
can also say that, it is the discontinuity of momentum of 
the ground state that leads to the jump of concurrence 
at the critical point. 

In order to investigate the role that the momentum 
plays on the change of various types of concurrence C<- a \ 
we reconstruct two degenerate ground states at the crit- 
ical point A as 

= cos ~hk> + e<* siii- life,), 
l* 2 ) = e-*+ sin 6 -\i>x)- cos t\fr), (23) 

where 9 6 [0,7r] and <j) £ [0, 2n]. For AT-site systems 
(N > 6), the pairwise concurrences of type a (total 
concurrence) of the two degenerate ground states |*i) 
and 1*2) are denoted as c[ Q ' and Cg (Cn and Ct%), 
respectively. A straightforward calculation shows that 
cj Q] = C [ 2 a] = Q for a > 1, while cf ] and C 2 1] are 
nonzero and depend on the parameters 9 and <f>. Ob- 
viously, states l^i) and 1*2) are no more the eigenstates 
of momentum in the general case. 

What we concern is the role the momentum plays on 
the jump of the concurrence. For N > 6 the difference 
between cj 1 ' and cJ, 1 ' is 



C [1 > ~ C l 2 ] 



£11 + £12 — £21 — £22! 



where 



-H(3Gij -1)(3Q„-1) 



Qij = 1 - 2xw + (-T^nXn cos ' 
+ (— y +:j XN cos sin 0, 



XN 



and 



H(x) 



N/2-2 



(s > 0) 

(£B < 0) 



(24) 



(25) 



(26) 



(27) 



(28) 



is the Heaviside step function. In the following, we will 
show that this difference reaches maxima at 9 — 0, ir for 
any <f> G [0, 2ir\. 

Notice that 3Qij — 1 > always holds when <f> = 
7r/2, 37r/2. Then Eq. (J2U can be rewritten as 



\C\ 



[i] 



C [ 2 ] \ = 3]&vxlr cos( 



(29) 



It is a monotonic function, which reaches its maxima at 

= O,7T. 



Now we prove that for any 9, the inequality 

IC^ -C [ 2 1] \<3\t NXN cos9\ (30) 

holds in the all range of (f>. It is convenient to consider 
the problem in the range 6,<p € [0, tt/2] without loss of 
generality. Since the functions sin 9, cos 9 and cos0 are 
positive in this range, we can prove the above conclusion 
in the following three cases: 
Case 1: 



3Qi2 
3Q 2 i 



Case 2: 



3Qi2 
3Q 2 i 



and 
Case 3: 



3Qi2 
3Q 2 i 



1 = A — 3xat cos </> sin > 0, 
1 = B — 3xjv cos sin 6> > 0; 



1 = A — 3\n cos sin < 0, 
1 = B — 3xArcos0sin0 > 0; 



1 = A — 3xn cos sin < 0, 
1 = B — 3\n coscj)sm9 < 0; 



where we have defined 

A = 2 - 6 X 2 N -'^nxI cos 9 >0, 
B = 2 - 6x1 + 3£nXn cos 9 >0, 
B > A. 

For all the above three cases we have 

3Qn — 1 = A + 3x./v cos </> sin # > 0, 
3Q22 — 1 = B + %xn cos ^ sin 6* > 0. 

Actually, in case 1, we have 

|cf-Cf| =3|6vx^cos0|. 
In case 2, from Ea. (|32() we have 

\D\ < |-4£ w x^cos0] , 



where 



D 



+XArCOS(/)sin6' 4- 2x N - 
Then the concurrence difference is 

ICP-Cfl = |£ll-£21-£ 22 | 

= - a \D\<2,\^ nXn cos9\ 
Similarly, for the case 3, we have 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



cj 1] -c ; 



Ml 



\£nXn cos9 \ < 3 kNX 2 N cos9\ . (39) 
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FIG. 4: The NN concurrence difference of the two recon- 
structed degenerate ground states for N = 6, 8, 12 and 16. 
Since all the rest types of concurrences are 0, the behavior of 
ICj 1 ' — CJj 1 '! and \Cti — Ct2\ are the same. It shows that the 
total concurrence difference reach the maxima when 9 — or 
7r for all <j> an d N . But the maxima decay exponentially with 
the size N of the system. 



For the cases that 9 and <j> are taken in the rest ranges, 
similar proof as presented in l|31 |) — l|39 [l can get the same 
conclusion. So from Eqs. and l|3UI) we conclude that 



the concurrence difference between C± and CJ, 11 reaches 
the maxima when 9 = 0, n for any tj>. 

In fact, for the cases N > 8 within the special domain 



of 9 ~ 0, 7r and any 
Hence we have 



it always hold that 3Qi 



1 > 0. 



|cf -C. 



Ml - 



(0-7T) 



(40) 



which reaches its maxima S^nXn a t ^ = Oj 71 "- Here 
the two combined states l^i) and |^) are just the 
eigenstates of momentum, i.e. the ground states at 
J = J c ±0+. 

As illustration, in Fig. 4 the concurrence differences, 
IqM _ as the functions of the parameters 9 and <j> 
are plotted for N = 6, 8, 12 and 16 systems. 

It shows that there always exists a maximal concur- 
rence difference, or a maximal concurrence jump, at the 
points where the two degenerate ground states are the 
eigenstates of momentum with k — 0, tt. Furthermore, in 
the thermodynamic limit, we have limAr^oo^jv = and 
limjv^ooX-/v = \i and then the maxima decays to zero 
exponentially The results indicate that the change of the 
momentum induces the maximal jump of |CW| or |CV| 
at the point A. So far, we can not say which types of 
concurrence is preferable to characterize the QPT, from 



FIG. 5: All types of concurrence differences of the degenerate 
ground states for N = 8 (a) and 10, 12 (b). Notice that the 
anomalous behavior in N — 8 case shows that only the total 
concurrence difference reach the maxima when 8 = or n for 
all (f> and N . 



the following arguments, we will see that an anomalous 
result for point B will give a final selection. 

Now we turn to consider the situation at the criti- 
cal point B. In this case, the differences of concurrence 
| _ £>H | anc [ \C T1 — Ct2 I can not be obtained analyti- 
cally. The numerical method is employed to calculate the 
differences of various types of concurrences. All types of 
concurrence difference as the functions of the parameters 
9 and <fi are plotted in Fig. 5. 

In Fig. 5(a), it shows that the pairwise concurrences 
of NN and NNN for N — 8 have jumps, while the rest 
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have no jumps. This result is different from that for the 
A point, in which case there is no jump for NNN concur- 
rence. Another interesting result is that the difference 
| cj 1 ' — C^ 1 ' | does not reach the maxima when the cor- 
responding two degenerate states are the eigenstates of 
momentum. This anomalous phenomenon indicates that 
the NN concurrence seems not to be sufficient charac- 
terize the QPT. However, Fig. 5(a) also shows that the 
difference \Cj-i — CV2I still obeys the same rule we ob- 
tained at point A. In Fig. 5(b) the corresponding results 
for N — 10 and 12 are plotted, which are similar to that 
of point A. Thus all the results imply that the differ- 
ence of total concurrence reach the maxima when two 
degenerate states are the eigenstates of momentum. In 
other words, it is the change of momentum that induces 
the jump of total concurrence. Based on all the results, 
we conclude that the total concurrence is a good candi- 
date to characterize the quantum critical behavior for the 
frustrated spin ring systems we concerned. 

V. DISCUSSION 

Summing up, in this paper we have shown how to es- 
tablish the connection between the concurrence jump and 
the change of groundstate momentum at the QPT critical 
point. All types of pairwise concurrence are investigated 
analytically and numerically. The results for both criti- 
cal points A and B indicate that the difference of total 
concurrence reach the maxima when the two degener- 
ate groundstates are just the eigenstates of momentum. 
It also reveal another interesting relation between cor- 
relation function and the concurrence. As mentioned in 
section II, the NN and NNN correlation function must 
have a jump at the energy-level crossing points, while 
the NN and NNN concurrences may not. But when the 
total concurrence is considered, it must have a jump at 
the critical points. However, the total correlation func- 
tion E^S, • Sj) g = i((S 2 - £. S?)) 9 (where S = £. S, 
is the total spin) has no jump at the critical points since 
the spins of the ground states are zero. 

Concerning the model for the case J\ 7^ J2 as illus- 
trated in Fig. 1, a straight forward calculation shows that 
the states j^i) and \4>2) defined in (|16|l are still the de- 
generate groundstates of the Hamiltonian if J1 + J2 = Jo- 
This means that J1 + J2 = Jo is the boundary of different 
quantum phases in the J\ — Ji plane. Starting from this 
observation we can extend our study to the more general 



case of J\ 7^ J2 to verify the conclusion obtained in this 
paper. It will appear in a successive paper. 
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